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ABSTRACT 

We study the comoving space density of X-ray-selected luminous active galactic nuclei (AGNs) and the 
obscured AGN fraction at high redshifts (3 < z < 5) in the Subam/XMM-Newton Deep Survey (SXDS) field. 
From an X-ray source catalog with high completeness of optical identification thanks to deep optical images, 
we select a sample of 30 AGNs at z > 3 with intrinsic (de-absorbed and rest-frame 2-10 keV) luminosities of 
Lx = 10 44-45 erg s" 1 detected in the 0.5-2 keV band, consisting of 20 and 10 objects with spectroscopic and 
photometric redshifts, respectively. Utilizing the 1 /V ma x method, we confirm that the comoving space density 
of luminous AGNs decreases with redshift above z > 3. When combined with the C/iant/ra-COSMOS result 
of Civano et al. (201 1), the density decline of AGNs with Lx = 10 44-45 erg s" 1 is well represented by a power 
law of (1 +z)~ 6 2±0 9 . We also determine the fraction of X-ray obscured AGNs with A^h > 10 22 cm -2 in the 
Compton-thin population to be 0.54^Q[g, by carefully taking into account observational biases including the 
effects of photon statistics for each source. This result is consistent with an independent determination of the 
type-2 AGN fraction based on optical properties, for which the fraction is found to be 0.59±0.09. Comparing 
our result with that obtained in the local Universe, we conclude that the obscured fraction of luminous AGNs 
increases significantly from z = 0toz>3bya factor of 2.5±1.1. 

Subject headings: galaxies: evolution — galaxies: active — X-rays: galaxies — surveys 



1. INTRODUCTION 

The evolution of active galactic nuclei (AGNs) carries key 
information for understanding the growth history of super- 
massive black holes (SMBHs) in galactic centers. Recent 
X-ray surveys have revealed that luminous AGNs have their 
number density peak at higher redshifts, z ~ 2-3, than less lu- 
minous ones (Ueda et al. 2003; Barger et al. 2005; La Franca 
et al. 2005; Hasinger et al. 2005; Silverman et al. 2008; Ebrero 
et al. 2009; Yencho et al. 2009; Aird et al. 2010). This be- 
havior is called "down-sizing" or "anti-hierarchical" evolu- 
tion. These results suggest that present-day massive SMBHs 
formed in the early epoch of the Universe, compared with 
less massive ones. A similar evolution is also observed in 
the star forming history of galaxies (e.g., Cowie et al. 1996; 
Kodama et al. 2004). This provides a further evidence for the 
cosmological "co-evolution" scenario of SMBHs and galax- 
ies, which is expected from the tight correlation between a 
SMBH mass and galactic bulge properties observed in the 
local Universe (see e.g., Magorrian et al. 1998; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). 

To elucidate the formation processes of SMBHs, it is very 
important to reveal the AGN evolution at the high redshifts be- 
fore the density peak, where the rapid growth of SMBHs took 
place. Optical/infrared AGN surveys, such as the Sloan Digi- 
tal Sky Survey, the Canada-France High-z Quasar Survey, and 
the UKIDSS Large Area Survey, are able to detect very lumi- 
nous quasars at high redshifts, now up to z — 7 (Jiang et al. 
2009; Willott et al. 2010; Mortlocket al. 2011), revealing that 
their space number density has a peak around z — 2 - 3 and 
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shows a rapid decline by 96% from z = 3 to z = 6 (Richards 
et al. 2006). Recently, the space density of lower luminos- 
ity AGNs has been also investigated from deep optical survey 
fields (Glikman et al. 201 1 ; Ikeda et al. 201 1), although the re- 
sults are still controversial due to the complexity of the correc- 
tions required for incompleteness and contamination. While 
these optical surveys are useful in detecting unobscured (so- 
called type-1) AGNs, they inevitably miss the population of 
obscured (type-2) AGNs, the major class of the AGN popula- 
tion (see e.g., Gilli et al. 2007). 

Hard X-ray observations at rest-frame energies above a few 
keV provide an ideal opportunity to detect obscured AGNs, 
thanks to their much smaller biases against absorption than 
the optical band and low contamination from stars in the host 
galaxies. Silverman et al. (2008) reports that the space den- 
sity of X-ray-selected luminous AGNs significantly declines 
at z > 3, similar to that of optically-selected QSOs, although 
they have to assume a completeness correction because the 
spectroscopic identification rate of the sample is not high 
(~ 50%). More recently, Brusa et al. (2009) and Civano et al. 
(2011) have studied X-ray detected AGNs at z ^ 3 from the 
COSMOS survey (Scoville et al. 2007; Hasinger et al. 2007; 
Elvis et al. 2009), one of wide and deep multi-wavelength sur- 
veys over a continuous area. A decline of the space density 
above z ~ 3 is confirmed for luminous AGNs with intrinsic 2- 
10 keV luminosities of >, 10 44 erg s . While the functional 
shape of this decline is found to be consistent with that of the 
exponential decline parameterization proposed by Schmidt et 
al. (1995), which is based on the optically-selected QSOs, it 
still has a large uncertainty due to the limited sample size. 
For quantitative comparison with the results obtained in other 
wavelengths and theoretical models, it is important to estab- 
lish the shape of the decline from a larger X-ray-selected AGN 
sample with high completeness. 

Another key observational property of AGNs is the type-2 
AGN fraction (or more precisely, the distribution of absorp- 
tion column density), which gives strong constraints on the 
environments around SMBHs including the obscuring tori and 
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host galaxies. According to the unified model of AGNs (An- 
tonucci 1993; Urry & Padovani 1995), the classification of 
AGNs depends only on the viewing angle between the ob- 
server and the rotation axis of the torus. It is well known that 
the type-2 AGN fraction decreases with increasing X-ray lu- 
minosity (Ueda et al. 2003; Steffen et al. 2003; La Franca et 
al. 2005; Hasinger 2008; Ueda et al. 201 1), indicating that the 
simplest unified model with a fixed geometry of tori does not 
hold. Several authors (e.g., La Franca et al. 2005; Treister 
& Urry 2006; Hasinger 2008) have furthermore reported that 
the type-2 fraction increases with redshift, suggesting an evo- 
lution of the structures around SMBHs. However, due to diffi- 
culties of type classification using photon-limited X-ray data 
and complex biases against obscuration, the intrinsic type-2 
AGN fraction at z > 3 has been only poorly investigated. 

In this paper, we investigate the evolution of the comoving 
space density of X-ray-selected luminous AGNs and type-2 
fraction at high redshifts (3 < z < 5) using the X-ray cata- 
log of Subaru/XMM-A/evvfon Deep Survey (SXDS), another 
wide and deep multi-wavelength survey. The organization of 
the paper is as follows. Section 2 describes the high-redshift 
AGN sample used in this study. In section 3, we explain the 
AGN type classification methods using X-ray or optical data. 
The results of the comoving space density and type-2 fraction 
are presented in section 4. Section 5 summarizes the conclu- 
sions. Throughout this paper, the cosmological parameters 
(H , ft M , ^a) = (70 km s" 1 Mpc" 1 , 0.3, 0.7) are adopted. The 
"log" symbol represents the base- 10 logarithm. Quoted errors 
denote those at la level. 

2. SAMPLE 

2.1. Identification of SXDS AGNs 

The Subaru/XMM-A/evvfon Deep Survey (SXDS; Furusawa 
et al. 2008; Ueda et al. 2008) is one of wide and deep multi- 
wavelength survey projects covering the radio to X-ray bands. 
The survey has an unprecedented combination of depth and 
sky area over a contiguous region of >1 deg 2 , centered at 
R.A. = 02 / T8 m and Dec. = -05''00'" (J2000). The main aims 
of the SXDS survey are to explore the nature of extragalactic 
populations over the whole history of the Universe, the co- 
evolution of galaxies and SMBHs, and the evolution of the 
large-scale structure, without being affected by cosmic vari- 
ance. 

XMM-Newton observed the SXDS field in seven pointings, 
with one deep (nominal exposure of 100 ks) observation in the 
center surrounded by six shallower (50 ks each) ones, in the 
0.2-10 keV band with the European Photon Imaging Cam- 
era (EPIC; Struder et al. 2001; Turner et al. 2001). From the 
combined images taken by three EPIC cameras (pn, MOS 1 , 
and MOS2) with a total area of 1.14 deg 2 , Ueda et al. (2008) 
detected 866, 1 1 14, 645, and 136 X-ray sources with sensitiv- 
ity limits of 6 x 10" 16 , 8 x 10" 16 , 3 x 10" 15 , and 5 x 10" 15 erg 
cm" 2 s" 1 in the 0.5-2, 0.5-4.5, 2-10, and 4.5-10 keV bands, 
respectively, and with a detection likelihood > 7 (correspond- 
ing to a confidence level of 99.91%). 

Akiyama et al. (2012, in preparation) summarize the results 
of the identification of the SXDS X-ray sources detected in 
the 0.5-2 keV and/or 2-10 keV band within the five com- 
bined fields-of-view of the Subaru/Suprime-Cam imaging, us- 
ing optical and infrared spectroscopic data as well as the 
deep imaging datasets obtained with UKIRTAVFCAM and 
Spitzer/IRAC. The 3a detection limit of the deep /?-band im- 
age reaches 27.5 magnitude. Thanks to the deep optical im- 
age, objects with X-ray to optical flux ratio, log(/x//R), up to 



+1 .3 can be detected, even for X-ray sources at the deepest X- 
ray detection limit in the soft X-ray band. The log(/x//R) dis- 
tribution of typical AGNs (-1 < log(/ x //R) < +1; Akiyama 
et al. 2000) is well enclosed within the detection limit. Fur- 
thermore, the 3a limit of the B-band image (28.5 magnitude) 
is deep enough to detect the sharp decline below Lya of high- 
redshift objects with red B-R color (>1.0 magnitude) even at 
the 7?-band detection limit. 

Among the 866 sources detected in the 0.5-2 keV from the 
whole SXDS field, 781 sources are also covered with the deep 
Subaru/Suprime-Cam optical images in B-, V-, R-, i-, and z- 
bands (Furusawa et al. 2008). In this study, we treat these 
781 sources located within the overlapping region of the X- 
ray and deep optical imaging data as a parent sample. Once 
candidates of Galactic stars, stellar objects that are close to 
a bright galaxy, and clusters of galaxies are removed, 733 
out of the 781 sources are remained as candidates of AGNs. 
Among them, 586 are spectroscopically identified with inten- 
sive optical and near-infrared spectroscopic observations in 
the field. For the other AGN candidates without spectroscopic 
identification, their photometric redshifts are estimated on the 
basis of photometric data in 15 bands covering far-UV to 
mid-infrared wavelength range. Akiyama et al. (2012) apply 
HyperZ photometric redshift code (Bolzonella et al. 2000), 
which determines a photometric redshift via template fitting 
with x 2 minimization, using both of galaxy and QSO Spec- 
tral Energy Distribution (SED) templates. 

In order to reduce the number of AGNs with photomet- 
ric redshifts significantly different from spectroscopic ones 
("outliers"), we apply two additional constraints in addition 
to the x 2 minimization, considering the properties of the 
spectroscopically-identified AGNs. First one is that the ob- 
jects with stellar morphology in the deep optical images are 
z > 1 broad-line AGNs. Almost all X-ray sources with stel- 
lar morphology are identified with type-1 AGNs at z > 1 in 
SXDS. They show bright nuclei and their observed optical 
lights are dominated by the nuclear components. This crite- 
rion may not be applicable to a survey with a deeper X-ray 
flux limit. Second one is the absolute magnitude range of the 
galaxy and QSO templates. Considering the absolute magni- 
tude range of spectroscopically-identified type-1 and type-2 
AGNs, we limit the z-band absolute magnitude ranges of the 
galaxy and QSO templates between -25.0 <M~< -20.0 and 
-26.5 < M z < -22.0, respectively. Further details of the pho- 
tometric redshift determination are discussed in Akiyama et 
al. (2012). 

The accuracy of the photometric redshifts is examined by 
comparing them with the spectroscopic ones. The median 
of ., A: , is 0.06 for the entire sample. We further exam- 

(l+Zspeo) r 

ine the accuracy by the normalized median absolute deviation 
(NMAD; a t ) following Brammer et al. (2008). For the en- 
tire sample, a z is 0.104, which is larger than that of the pho- 
tometric redshift estimations for X-ray-selected AGNs with 
medium band filters (Cardamone et al. 2010; Luo et al. 2010). 
The a z for type-1 AGNs (0.201) is larger than that for type- 
2 AGNs (0.095). This is because there is no strong feature 
in the SEDs of the type-1 AGNs except for the break below 
Lya. However, since the break can be covered by the deep 
optical images for AGNs at z > 3, the discrepancy between 
z S pec and Zphot is smaller for AGNs at z > 3. Regarding the 
completeness and contamination for the z > 3 AGN sample, 
one AGN at z spe c < 2.9 shows z p hot > 3.0. In contrast, five 
AGNs at z spec < 2.9 have z phot > 3.0. 
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FIG. 1 . — Distributions of SXDS AGNs on the R — z and B—R color-color diagram (left: z > 3 objects, right: z < 3 objects). The red squares and green circles 
represent type- 1 and type-2 AGNs, respectively. AGNs with spectroscopic and photometric redshifts are plotted with filled and open symbols, respectively. The 
arrows illustrate objects whose B—R colors are their lower limits. The dashed line displays a possible criterion for z > 3 AGNs with dropout selection. 



It needs to be emphasized that in order to determine photo- 
metric redshifts of z > 3 AGNs reliably, deep optical imaging 
data are crucial since they can detect sharp decline of flux be- 
low the wavelength of rest-frame Lya emission. Finally, most 
of the X-ray sources are identified as AGNs at well estimated 
redshifts with a very high degree of completeness (99.2%). 
Photometric redshifts cannot be determined for six objects 
that are detected in less than 4 bands. We discuss their possi- 
ble contributions in sections 2.3 and 4.1. 

Dropout method is a simple selection criterion for galax- 
ies and AGNs at high redshifts, and sometimes such selection 
is preferred than the selection with the photometric redshifts, 
because the photometric redshifts can be affected by a catas- 
trophic failure (Aird et al. 2010). In order to compare the 
dropout selection with the photometric redshift selection, we 
plot R - z and B-R color-color diagram for AGNs at z > 3 
(left) and z < 3 (right) in figure 1. The z > 3 AGNs show 
red B-R color, which is consistent with the break below Lya 
emission. The red color can be a crude check of the photomet- 
ric redshifts, which are determined by the SEDs with all 15 
bands. Most of the z > 3 type-1 AGNs have blue R-z color as 
well, and they are enclosed within the criterion (dotted lines) 
that is determined to exclude most of the z < 3 AGNs (three 
z < 3 type-2 AGNs that are enclosed within the criterion have 
photometric redshifts of 2.6-2.8, close to z = 3). In contrast, 
z > 3 type-2 AGNs have redder R-z color than type- 1 AGNs 
on average, and some of them have similar colors to AGNs 
at lower redshifts. Thus, selection based on the optical colors 
only is useful in general but may miss a part of type-2 AGNs 
that have red UV color above Lya emission. 

2.2. z > 3 AGN Sample 

For our study, we select a high-redshift AGN sample at 
z > 3 detected in the 0.5-2 keV band, on the basis of spectro- 
scopic and photometric redshifts determined by Akiyama et 
al. (2012). Thanks to the K-correction effect, even absorbed 
AGNs can be observed in the 0.5-2 keV band (corresponding 
to the rest-frame energies above 2-8 keV for z > 3), where the 
maximum sensitivity of XMM-Newton is achieved. Our final 
z > 3 AGN sample consists of 30 objects, among which 20 
and 10 AGNs have spectroscopic and photometric redshifts, 
respectively. The properties of the sample are summarized in 



table 1 and the optical spectra of the spectroscopic sample are 
displayed in figure 2. Figure 3 exhibits the distribution of red- 
shift for the z > 3 AGN sample. Possible effects related to the 
uncertainties in the photometric redshifts are discussed in the 
next subsection. 

2.3. Uncertainties of Photometric Redshifts 

Since the photometric redshifts are determined only by 
broad-band photometries, they could have uncertainties un- 
like spectroscopic redshifts. Figure 4 illustrates examples of 
the normalized probability distribution functions (PDFs, P(z)) 
for the photometric redshift for two sources without spectro- 
scopic redshifts, SXDF0154 and SXDF0321, calculated from 
the distribution of the \ 2 value obtained in the SED fitting. 
SXDF0154 has a primary solution (i.e., the highest peak of 
PDF) at z > 3, and hence is included in our z > 3 AGN sam- 
ple. By contrast, SXDF0321 is not included in the sample 
since the PDF has its peak at z < 3. Both of them, however, 
have finite probabilities at z < 3 and z > 3. 

To take into account the shape of the PDF in each object 
statistically, we consider the contribution from all 0.5-2 keV 
selected AGNs in the Akiyama et al. (2012) catalog to the 
z > 3 AGN sample, using the normalized PDFs between z = 
and z = 5. For the sources with spectroscopic redshifts, we 
define their PDFs as delta-functions peaked at their measured 
redshifts. It is found that the "effective" number of objects in 
the z > 3 AGN sample, computed by summing the fraction of 

the normalized PDFs in the z = 3-5 region (J^P(z)dz) from all 
objects, is 30.6, when the six unidentified objects are ignored. 
This value is almost the same as the number of the best-fit 
z > 3 AGN sample described in the previous subsection. Even 
if we include the unidentified objects for the calculation, the 
effective number is computed to be 33.7, which is still within 
the la error of the number of our z > 3 AGN sample. We also 
confirm that there are no significant changes in the luminos- 
ity and absorption column density distributions of the sample 
when we consider the PDFs of photometric redshifts; here the 
contribution of each object to the distributions is added with a 
weight of P{z) over the redshift range of z = 3-5. Hence, we 
hereafter consider only the best-fit values of the photometric 
redshifts. 
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FIG. 2. — Optical spectra of the twenty z > 3 AGNs with spectroscopic redshifts. The vertical axis represents the flux density per wavelength (some spectra 
are not flux-calibrated). The vertical dashed lines in each panel indicate the positions of emission lines typically seen in AGNs. The gray hatched regions exhibit 
wavelength ranges affected by strong night sky lines or atmospheric absorption lines. 



3. ANALYSIS 

3.1. Classification of AGNs 

For our discussion of the obscured AGN fraction, we need 
to classify each object in our sample as a type-1 or type-2 
AGN. There are two different ways to classify AGN: one is 
based on the optical spectra (optical type) and the other is 
based on the absorption column densities, A^h, in the X-ray 
spectra (X-ray type). The optical type is basically determined 
by the presence of broad, permitted emission lines. Con- 



sidering the distribution of the full width at half maximum 
(FWHM) of the broad-line AGNs in the local universe (Stern 
& Laor 2012), we classify our spectroscopic sample: those 
showing significant broad lines with velocity widths larger 
than 1,000 km s" 1 are optical type-1 and the others are op- 
tical type-2. In the z > 3 sample, all but one (SXDF0809) 
broad-line AGNs have the FWHM of the broad-line width 
larger than 2,000 km s" 1 . Thus, their optical classifications 
are not changed even if we use the threshold of 2,000 km s" 1 
instead of 1,000 km s" 1 . The CIV line of SXDF0809 has the 
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FIG. 2. — Continued. 



FWHM of 1,470 km s _1 . Although it is classified as a type-1 
AGN since the CIV line shows broader wing, it needs to be 
noted that it also shows a strong Hell line, which is typical 
for narrow-line AGNs. The contribution from the narrow-line 
component is large in this object. 

For the 10 AGNs without spectroscopic redshifts (i.e., those 
that do not show any significant line features or do not have 
optical/near-infrared spectra available), we classify them by 
the type of a template spectrum (i.e., QSO or galaxy). For 
spectroscopically-identified AGNs, there is a good correla- 



tion between the spectral type and the best-fit SED type: 
type-1 (type-2) AGNs are well fitted with QSO (galaxy) tem- 
plates. Specifically, among the twenty z > 3 AGNs with 
spectroscopic identification, the spectroscopic and photomet- 
ric types match well for 15 objects. We therefore classify 
spectroscopically-unidentified objects fitted better with QSO 
(galaxy) template into type-1 (type-2) AGNs. Considering 
the mismatch rate of 25%, we estimate that less than two (= 
9*0.25 - 1*0.25) out of the nine AGNs with galaxy spectra 
could be truly type-1 AGNs; this uncertainty will be taken 
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into account in the discussion on optical type-2 AGN fraction 
(section 4.2.2). We do not consider combined templates with 
galaxy and AGN in the photometric redshift determination. It 
will not affect the classification of the z > 3 AGNs, since at 
z > 3 the SXDS sample covers only luminous AGNs and the 
SEDs of luminous type-1 AGNs are completely dominated by 
QSO components. 

A general problem of relying on the optical type is that the 
classification may depend on the quality of the available op- 
tical spectra, since good signal-to-noise ratio is required to 
detect less-luminous broad-emission lines on top of stellar 
continuum emission of host galaxies. Nevertheless, such ef- 
fects are not expected to be significant in our sample since it 
consists of luminous AGNs for which contamination from the 
host galaxies is negligible. 

X-ray absorption is an alternative good indicator of AGN 
type, as expected from the unified scheme (Antonucci 1993; 
Urry & Padovani 1995) : type-2 AGNs are viewed through the 
dust torus, and hence have higher absorption column densities 
than type-1 AGNs. In the present work, we adopt A^h = 10 22 
cm" 2 as the dividing criterion: AGNs with log A^h < 22 and 
> 22 are classified as X-ray type-1 and type-2, respectively. 
This criterion is adopted by many authors, and is known to be 
generally in good agreement with the optical type (see e.g., 
Ueda et al. 2003). The merit of using X-ray classification is 
that observational biases in X-ray surveys can be well defined 
and can be corrected in a straightforward way. Thus, we basi- 
cally adopt the X-ray type for our study, although consistency 
with results based on the optical type is also discussed. 

3.2. Absorption Column Density 

To derive the absorption column density (A^h) for each 
AGN, we use an X-ray hardness ratio (HR) between two en- 
ergy bands, since it is difficult to perform spectral fitting due 
to the limited photon statistics for most of the X-ray sources 
in our sample. We define the HR value as 



UR = (H-S)/(H+S) 1 



(1) 



where S and H are the vignetting-corrected count rates 5 in 
the 0.3-1.0 keV and 1.0-4.5 keV bands, respectively. We find 
that this choice of bands gives the best sensitivity for the de- 
termination of the A^h value for A^h ~ 10 22 cm -2 , for objects 
at z = 3-5. Since these count rates are not available in the 
Ueda et al. (2008) catalog, we re-analyzed the XMM-Newton 
images in these two bands at the fixed source positions, utiliz- 
ing exactly the same datasets and analysis software as used in 
Ueda et al. (2008). The relationship between HR and redshift 
of our AGN sample is shown in figure 5. 

Figure 6 shows the relationship between HR and column 
density at z = 3, 4, and 5. Here and after, we assume a power 
law with a photon index of 1 .9 plus a reflection component 
from cold matter with a solid angle of 2ir as the intrinsic spec- 
trum, following Ueda et al. (2003). We derive the column 
densities of our AGN sample from the redshift and observed 
HR value using this relation, and then calculate their intrinsic 
(de-absorbed and rest-frame 2-10 keV) fluxes and luminosi- 
ties. Where an observed HR corresponds to log A^h < 20, we 
set it to log A^h = 20. The column density, flux, and luminosity 
are also summarized in table 1, and the flux distribution is pre- 
sented in figure 7. Figure 8 displays the relationship between 

5 corresponding to those of the pn detector, obtained by dividing the com- 
bined pn+MOSl+MOS2 image by the pn-equivalent exposure map; for de- 
tails see Ueda et al. (2008) 
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FIG. 5. — Hardness ratio (HR) versus redshift plot for our AGN sample. The 
red and green symbols correspond to type-1 and type-2 AGNs, respectively. 
The filled and open circles represent spectroscopic and photometric redshifts 
samples, respectively. The magenta dashed line corresponds to log N\\ = 22. 
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FIG. 6. — The relationship between the hardness ratio HR and absorption 
column density N^. A power law with a photon index of 1.9 plus a reflection 
component from cold matter with the solid angle Q = 2-7T is assumed as the 
intrinsic spectrum. The red solid, green dashed, and blue dot-dashed lines 
correspond to the redshifts of 3, 4, and 5, respectively. 



the intrinsic luminosity and redshift for our AGN sample. The 
magenta curves represent the detection limits for AGNs with 
log A^h = 20, 22, and 23, calculated from the 0.5-2 keV count 
rate at which the survey area becomes one tenth of the maxi- 
mum value. As noticed from the figure, the survey is incom- 
plete for sources with log Nn > 23 in the lowest luminosity 
range below log Lx ~ 44.2 at z > 3.5. However, we con- 
firm that this does not affect our estimate on the space density 
of all Compton-thin AGNs with log A^h < 24 over the statis- 
tical errors as far as the fraction of heavily obscured AGNs 
with A^h = 23-24 is not extremely high. We find that all of 
our AGN sample have intrinsic luminosities of log Lx > 44 
erg s -1 , and hence belong to the QSO class. While the X-ray 
AGN types match the optical ones for 18 out of 30 AGNs in 
our sample, they do not for the others. This is largely due to 
the statistical error of column density. We can, however, cor- 
rect for it by considering a conversion matrix from an intrinsic 
to observed A^h distributions (see section 4.2.1). 

4. RESULT AND DISCUSSION 
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FIG. 7. — Flux distribution for z > 3 AGN sample. The definition of the 
symbols is the same as in figure 3. 
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FIG. 8. — Correlation between the intrinsic (de-absorbed and rest-frame 
2-10 keV) luminosity and redshift for our AGN sample. The definition of the 
symbols is the same as in figure 5. The magenta lines represent the detection 
limits for AGNs with log Nu = 20, 22, and 23, from bottom to top, calculated 
from the 0.5-2 keV count rate corresponding to one tenth of the total survey 
area. 

4.1. Evolution of the Comoving Space Density of AGNs at 

z>3 

To investigate the cosmological evolution of AGNs at z > 3, 
we calculate the comoving space density from our sample uti- 
lizing the 1 /V max method (Schmidt 1968). This method takes 
into account the fact that a fainter source can be found in a 
smaller volume because of the survey area dependence on flux 
(or more correctly, count rate). The maximum available vol- 
ume V max is calculated by the following formula: 

X n(c(L x ,z,N H W +zfd A (z) 2 Az)dz, (2) 

dz 

where tt(c(Lx,z,Nn)) is the sky coverage at the count rate 
c(L x ,z,Nh) expected from a source at the redshift z with the 
intrinsic luminosity Lx and the absorption column density Nu, 
cIa and c^p are the angular distance and the look-back time 
distance per unit redshift, respectively, both of which are func- 
tions of z- Here z m j n and z max are the lower and the upper edges 
of the redshift bin used for computing V m3X , respectively. In 
the present paper, we separate the redshift interval of 3 < z < 5 
into five bins (z = 3.0-3.2, 3.2-3.4, 3.4-3.8, 3.8-4.3, and 4.3- 
5.0), and calculate V max for each AGN in the corresponding 
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FIG. 9. — Comoving space density of X-ray-selected, luminous AGNs 
plotted as a function of redshift. The red squares and blue triangles indicate 
our results at log Lx > 44 and those of Civano et al. (201 1) at log Lx > 44. 15, 
respectively. The errors shown are Icr. The gray dashed arrows represent the 
strict upper limits at each redshift bin (for visuality, they are slightly shifted 
in the x-axis direction), calculated by assuming that all of the unidentified 6 
objects are included in the bin. The magenta curve corresponds to the best-fit 
model fitted to the SXDS data with the form of (1 +z)~ 6 ' 2 , obtained from the 

re^ftluftlDin! 1 After calculating Vmax, we sum up the reciprocal 
value in each redshift bin: 



•= E 



i 



v m 



(3) 



where <f> is the comoving space density in the range z = Zmin - 
Zmax, and i is the index of the sample AGN falling into the 
redshift bin. The la error is estimated as: 



(70 = 



\ 



n<Z<Zn 

E 



i 



v m 



(4) 



Where only one source occurs in a given redshift bin, we 
adopt the Poisson upper and lower limits corresponding to 
the confidence range of Gaussian lcr formulated by Gehrels 
(1986). 

Figure 9 displays, as a function of redshift, the comoving 
space density of AGNs with luminosities log Lx = 44^45 erg 
s -1 obtained from the SXDS sample (red squares). The gray 
dashed arrows represent the strict upper limits at each redshift 
bin, calculated by assuming that all of the unidentified 6 ob- 
jects are included in the bin. The results from the COSMOS 
survey (Civano et al. 2011) for AGNs with log L x > 44.15 
erg s -1 are also plotted for comparison (blue triangles). As 
can be seen, our results are consistent with those of COSMOS 
within the statistical errors, indicating a significant decline in 
the AGN space density from z = 3 to higher redshifts. 

When we fit the SXDS result with a power-law form of 
oc (1 +z) p3 , we obtain p3 = -7.0 ± 1.8. We confirm that the 
index does not change over the error even when we consider 
the PDFs of photometric redshifts for all objects including 
the six unidentified ones (see section 2.3). A simultaneous 
fit to the SXDS and COSMOS data yields a best-fit slope of 
p3 = -6.2 ±0.9, which is shown by the solid magenta curve in 
figure 9. We find that an exponentially decaying model above 
a critical redshift z C 2 = 2.7, with the form 10 p4(z_Zc2) as adopted 
by Schmidt et al. (1995), also describes these data well; we 
obtain p4 = -0.65 ±0.15 (SXDS only) and p4 = -0.60 ± 0.08 
(SXDS+COSMOS). These results on p3 or p4 give the tight- 
est constraints on the decay profile of the space density of 



X-ray-selected luminous AGNs with log Lx = 44^45 erg s -1 
by utilizing a total of 80 objects. The parameters are little af- 
fected by filtering the SXDS sample with log L x > 44.15, the 
same criterion for the COSMOS sample. 

We confirm that the shape of the space-density evolution 
of X-ray-selected luminous AGNs is consistent with that de- 
rived from optical QSO surveys within current uncertainties. 
With the exponential model, Schmidt et al. (1995) obtained 
p4 = -0.43 ± 0.04 from optically-selected QSOs with M B < 
-26, which roughly corresponds to log Lx > 44.5 assuming 
a typical SED of QSOs (see e.g., Ueda et al. 2003). More 
recent results by Richards et al. (2006) give p4 ~ -0.67 for 
optical QSOs at a magnitude limit of M, < -27.6. Our result 
(p4 = -0.60 ± 0.08) obtained from X-ray surveys for AGNs 
with log Lx > 44 lies between the two optical results, as al- 
ready noted by Brusa et al. (2009). Since X-ray-selected sam- 
ples contain both type-1 and type-2 AGNs, the same decline 
profile as that of optically-selected (type-1) QSOs suggests no 
significant cosmological evolution of the type-2 AGN fraction 
above z > 3, at least in the bright luminosity range. This re- 
sult is consistent with the conclusion by Hasinger (2008) that 
the type-2 AGN fraction increases with z but saturates beyond 
z = 2 orz= 3. 

4.2. Fraction of Type-2 AGNs at z = 3-5 

To investigate the cosmological evolution of the fraction of 
type-2 AGNs (/type-2) in the whole Compton-thin population 
(i.e., with Nh < 10 24 cm" 2 ) from z = to z = 5, we constrain 
/type-2 at z = 3-5 in the intrinsic luminosity range of log Lx 
= 44-45 using our SXDS sample. We assume no redshift 
evolution within the redshift range of z = 3-5, following the 
previous discussion. There are observational biases in the ob- 
served ratio of type-1 and type-2 AGNs, which require care- 
ful correction to obtain the true fraction. As mentioned in 
section 3.1, we primarily adopt the "X-ray type" defined by 
the absorption seen in the X-ray spectra, for which an eval- 
uation of such biases can be made without assuming a rela- 
tionship between the optical and X-ray properties. We finally 
discuss the cosmological evolution of the type-2 AGN frac- 
tion by comparing our work with previous results obtained at 
lower redshifts. 

4.2.1. X-ray Type-2 Fraction 

In our analysis, we consider the two effects. (1) For a 
given intrinsic luminosity, X-ray type-2 AGNs become less 
detectable due to the reduction of the observed count rate by 
photoelectric absorption. (2) Due to the limited number of 
photons detected, there is a statistical fluctuation of the best- 
fit column density derived from the observed hardness ratio 
(i.e., an AGN can apparently have a different absorption col- 
umn density from the true value). It must be noted that the ob- 
served column densities are not simply subject to Gaussian or 
Poissonian errors but have complex probability distributions 
with unsymmetric negative and positive errors. This could 
make the "observed" distribution simply based on the best-fit 
column density of each object quite different from the "intrin- 
sic" one. 

To obtain the true absorption distribution by taking into ac- 
count these two biases simultaneously, we consider a "con- 
version matrix" from the intrinsic Nh distribution into the ob- 
served one, similar to that used by Ueda et al. (2003). Here we 
consider the X-ray luminosity function (XLF), the co-moving 
number density per unit co-moving volume per log Lx inter- 
val as a function of Lx and z, to predict the expected number 
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of AGNs with a given range of A^h from the SXDS survey, 
using the area curve given as a function of count rate (Ueda et 
al. 2008). Then, by simulating the statistical fluctuation in the 
observed count rates (and hence hardness ratios) according to 
the real survey characteristics, we reproduce the expected dis- 
tribution of the best-fit jVh from each input A^h value. Finally, 
we perform a fit to the observed histogram of A^h with the 
model prediction calculated via the matrix, by changing the 
type-2 fraction as a free parameter. With this procedure, we 
can most reliably constrain the "intrinsic" absorption distribu- 
tion in a statistical sense, even if it is hard to determine A^h of 
an individual object. 

For the shape of the 2-10 keV XLF, we adopt the 
luminosity-dependent density evolution (LDDE) model of 
Ueda et al. (2003) by further introducing a power-law decline 
of the comoving space density at high redshifts as reported in 
subsection 4.1. Thus, the evolution factor in equation (16) of 
Ueda et al. (2003) is replaced by the form: 

(l+z) pl 



e(z,L x ) = { 



l+z 



p2 



(1+Zc) 



l+Zc) 



1+Z 
i+Zc2 



p3 



[Z>Zc2~], 



where pi = 4.23, pi = -1.5, and z c = 1.9 are adopted. 



(5) 
fol- 



lowing Ueda et al. (2003). Here the cutoff redshift, beyond 
which the comoving space density of AGNs decreases, is 
set to z C 2 = 2.7, following Schmidt et al. (1995). We adopt 
p3 = -6.2 for the power-law slope at z > z c i based on our re- 
sult from both the SXDS and COSMOS samples. We also use 
the 'W H function" introduced by Ueda et al. (2003), f(N n ), i.e. 
the normalized probability distribution function for the ab- 
sorption column density. In our calculation, f(Nn) is assumed 
to be independent of luminosity and redshift for simplicity. 
We confirm that the effect caused by the luminosity depen- 
dence of the type-2 fraction (e.g., Ueda et al. 2003; Hasinger 
2008) on our result is smaller than the statistical uncertainty, 
because the luminosity range of our sample is narrow (log L x 
= 44-45). 

The expected number of AGNs detected, N, in the intrin- 
sic luminosity range log Lx = 44-45, z = 3-5, and intrinsic 
column density range A^h = A^H.min-A^H.max is computed as 

N= d\ogN H / c/logLx / dz 

f(N H )^^-n(c(L x ,z,N H ))(l+z) 3 d A (zfAz). 
dlogLx dz 

(6) 

The N values, calculated by assuming a flat A^h function, give 
the relative detection efficiency for AGNs with different col- 
umn densities. Figure 10 shows the result, normalized to unity 
at log A^h =20.0-20.5. 

Figure 1 1 illustrates the conversion matrix from the intrin- 
sic to observed A^h distributions, calculated for a flat /(Nh) in 
the range of log A^h = 20-24. Each panel of figure 1 1 repre- 
sents the expected A^h distribution for AGNs with input intrin- 
sic column densities of log A^h = 20-21, log A^h = 21-22, log 
A^h = 22-23, and log A^h = 23-24. Here, it is again empha- 
sized that the statistical errors in the count rates (eventually, 
those in the best-fit values) are taken into account, based 
on the observed correlation between the count rate and its er- 
ror, from the re-analysis described in section 3.2. It is seen 
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FIG. 10. — Detection efficiency as a function of intrinsic absorption column 
density, normalized to unity at log Nu = 20.0-20.5. 
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FIG. 1 1 . — Conversion matrix from the intrinsic 7Vh distribution into the 
observed one. The assumed intrinsic Nh distribution is plotted in the top 
right of each panel. 

from figure 1 1 that the total number of detected AGNs de- 
creases with the input A^h and the observed A^h distribution 
has a wide dispersion with respect to the input A^h value due 
to statistical fluctuations. 

Using this matrix, we perform a Poisson maximum likeli- 
hood fitting to the observed A^h histogram of our AGN sam- 
ple. We assume a flat shape for the A^h function in two A^h 
regions, log A^h = 20-22 and log A^h = 22-24, with the type-2 
AGN fraction as a free parameter to be determined through 
the fitting. The shape of the XLF is fixed to the model de- 
scribed above, while its normalization is tuned to reproduce 
the total number of observed AGNs at z = 3-5 for our best-fit 
A^h function. It is confirmed, however, that the choice of the 
XLF model does not affect our conclusion, by introducing the 
luminosity and density evolution (LADE) model (Aird et al. 
2010) instead of the LDDE model. We finally obtain a X-ray 
type-2 AGN fraction of 0.54^[J |g from our sample. Figure 12 
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FIG. 12. — Observed histogram of 7Vh with la Poisson errors (black 
squares). The red solid line represents the prediction from the /Vh function 
with the best-fit type-2 AGN fraction. 

shows the observed histogram of A^h with la Poisson errors 
(black squares), on which the model prediction from the 
function with the best-fit type-2 AGN fraction is overplotted 
(red solid line). 

4.2.2. Optical Type-2 Fraction 

To confirm the result obtained in the previous subsection, 
here we also estimate the type-2 AGN fraction based on the 
optical type of each AGN, not by referring to the X-ray ab- 
sorption. In our sample, we have a total of 17 optical type-1 
and 13 optical type-2 AGNs. To derive the intrinsic optical 
type-2 AGN fraction, we have to assume a relationship be- 
tween optical type and absorption column density, since the 
latter determines the observational biases in an X-ray survey. 
Here, we simply assume that our optical type-1 and type-2 
AGNs correspond to those with log A^h < 22 and log Nh = 
22-24, respectively, following Ueda et al. (2003). Then, we 
can utilize the detection efficiency plotted in figure 10 to cor- 
rect for the detection bias. Since the optical properties are in- 
dependent of the X-ray spectral information, we do not need 
to consider the statistical error in A^h in contrast with the pre- 
vious case. In this way, the optical type-2 fraction for AGNs 
with log Lx = 44^-5 erg s" 1 at z = 3-5 is estimated to be 
/type-2 = 0.59±0.09. Even when we consider the possible un- 
certainties in the classification of the 10 photometric redshift 
sample (section 3.1), the type-2 AGN fraction can be at least 
0.52±0.09. These values are consistent with the X-ray re- 
sult. The statistical error becomes smaller than that for the 
X-ray type-2 fraction because there is no coupling between 
the two types caused by the statistical fluctuation in the A^h 
value. We must note, however, that the correspondence of the 
optical type to X-ray absorption may have some ambiguities 
that could produce additional systematic uncertainties. 

4.2.3. Evolution of Type-2 Fraction in Luminous AGNs 

By fully taking into account the observational biases, in- 
cluding the photon statistics of the X-ray survey data, we have 
derived the X-ray type-2 AGN (log A^h > 22) fraction in the 
Compton-thin population to be 0.54^q'[9 at z = 3-5 and in the 
luminosity range of log Lx = 44-45. The consistency with the 
optical type-2 fraction of 0.59±0.09, which is independently 
measured, supports the robustness of our result. Recently, 
Masters et al. (2012) roughly estimate the type-2 AGN frac- 
tion with log Lx > 44. 15 at z ^3-5 to be ^75%, by comparing 



the normalizations of X-ray (type-1 + type-2) and UV (type-1 
only) luminosity functions derived from the COSMOS data. 
Their estimate is somewhat larger than our result, although it 
should be very sensitive to the adopted a ox value used for the 
luminosity conversion. 

We can now compare our X-ray type-2 AGN fraction with 
that obtained at lower redshifts to investigate its cosmologi- 
cal evolution. Recently, Burlon et al. (2011) have compiled 
an AGN sample in the local Universe detected in the 3-year 
Swift/BAT survey performed in the 15-55 keV band, where 
no detection bias is expected for Compton-thin AGNs (log 
A^h < 24). We find the X-ray type-2 fraction in their sample 
to be 0.22 ± 0.06 in the luminosity range as corresponding to 
log Lx = 44-45 in the 2-10 keV band, converted assuming a 
photon index of 2. Our result at z > 3 is significantly higher 
than the local value, at a ?»90% confidence level based on a 
X 2 test, by a factor of 2.5 ± 1.1 (the error is la). When we in- 
stead adopt the optical classification for the SXDS sample, the 
factor becomes 2.7 ±0.9 and the confidence level corresponds 
to 99.9%. Thus, we establish the evolution of type-2 fraction 
from z = to z > 3 in the high-luminosity range, where it has 
not been well constrained so far. 

We confirm that our results are consistent with the formu- 
lation of Hasinger (2008) for the type-2 fraction given as a 
function of luminosity and redshift. He finds that the type-2 
fraction normalized at log Lx = 43.75 (to correct for the lumi- 
nosity dependence) increases as (1 +z) 62 and then saturates 
above z — 2. The formula given in Hasinger (2008) gives 

/type-2 = (-0.11 ±0.12) * (log L x - 43 .75) + (0.578 ±0.081) 

(z= 3.2 -5.217) 

By substituting log Lx = 44.6 as the averaged luminosity in 
our sample, the type-2 fraction is calculated to be 0.49±0.13 
at z = 3.2-5.2. Note that Hasinger's "type-2" definition uti- 
lizes both optical and X-ray information, which is not exactly 
the same as ours; in fact, he argues that the criterion roughly 
corresponds to log A^h = 21.5, instead of 22, and therefore we 
have to reduce / yp e-2 to match our definition of X-ray type-2 
fraction. If we assume a flat A^h function above log Nu = 21.5, 
the value of type-2 AGN fraction in his definition is reduced 
by 20% in our definition; specifically, /t ype -2 = 0.49 ±0.13 
in his definition corresponds to an X-ray type-2 fraction of 
0.39±0.10 in ours. Thus, the formula of Hasinger (2008) is 
in agreement with the X-ray type-2 fraction derived from our 
SXDS sample within the errors, although our results may fa- 
vor a slightly higher type-2 fraction at z > 3. 

5. CONCLUSIONS 

Thanks to the deep and wide area coverage of the SXDS 
survey with extensive multi-wavelength follow-up observa- 
tions, we are able to construct a highly complete X-ray- 
selected AGN sample at high redshifts (3 < z < 5), consist- 
ing of 20 and 10 AGNs with spectroscopic and photometric 
redshifts, respectively. The conclusions are summarized as 
follows: 

• We derive the comoving space density of Compton- 
thin AGNs at z = 3-5 with intrinsic luminosities of Lx 
_ jq44-45 er g s -i at z > 3, using the 1/V ma x method. 
We confirm the trend that the space density signifi- 
cantly declines with redshift as reported in previous 
works. Combining our SXDS result with that from 
the COSMOS survey (Civano et al. 201 1), we find that 
the comoving space density decreases as (1 + z )- 6 - 2 ± - 9 
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or 10 -(0.60±0.08)Xfe-2.7) i providing so .f ar the best con . 

straint on the declining profile of X-ray-selected AGNs 
atz > 3. 



By carefully taking into account observational biases, 
including statistical uncertainties in the absorption col- 
umn densities derived from the X-ray data, we derive 

among 



the obscured AGN fraction with Nh > 10 22 cm 2 



those with N H < 10 24 cm" 2 to be 0.54^ at z = 3-5 



in the luminosity range of L x = 10 erg s . From 
comparison with the results obtained at lower redshifts, 
we establish that there is a significant redshift evolution 
of obscured fraction from z = 0toz>3bya factor of 



2.5 ± 1 . 1 in the high-luminosity range. 
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TABLE 1 

Properties of our high-redshift AGN sample. 
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a 0: photometric redshift AGNs, 1: optical type-1 AGNs, and 2: optical type-2 AGNs 

^SED models used for the calculations of photometric redshifts: 1 ; QSO template and 2; Galaxy template 

c Chosen redshift; spectroscopic redshift; photometric redshift 

d Hardness ratio defined as "HR = (H + S)/(H-Sy\ S: 0.3-1.0 keV count rate and H: 1.0-4.5 keV count rate 
e The absorption column density in units of cm -2 . 

^The intrinsic (de-absorbed and rest-frame 2-10 keV) flux in units of 10~ 15 erg cm -2 s -1 . 
g The intrinsic (de-absorbed and rest-frame 2-10 keV) luminosity in units of 10 44 erg s -1 . 



